GRB 090423 at a redshift of z ~ 8.1 
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Gamma-ray bursts (GRBSs) are produced by rare types of massive stellar explosions. 
Their rapidly fading afterglows are often bright enough at optical wavelengths, that they 
are detectable up to cosmological distances. Hirtheto, the highest known redshift for a 
GRB was z = 6.7 (ref. 1), for GRB 080913, and for a galaxy was z = 6.96 (ref. 2). Here we 
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report observations of GRB 090423 and the near-infrared spectroscopic measurement of 
its redshift z = 8.1I0 3. This burst happened when the Universe was only ~ 4% of its current 
age 3 . Its properties are similar to those of GRBs observed at low/intermediate redshifts, 
suggesting that the mechanisms and progenitors that gave rise to this burst about 600 
million years after the Big Bang are not markedly different from those producing GRBs 
~ 10 billion years later. 



GRB 090423 was detected by NASA's Swift satellite on 23 April 2009 at 07:55:19 UT as a double-peaked 
burst of duration T90 = 10.3 ±1.1 s. As observed by Swift's Burst Alert Telescope (BAT) 4 , it had a 
15-150 keV fluence F = (5.9 ± 0.4) x 10 erg cm and a peak energy E p = 481 6 , keV (errors at 90% 
confidence level). Its X-ray afterglow was identified by Swift's X-ray Telescope (XRT), which began 
observations 73 s after the BAT trigger 5 . A prominent flare was detected at t ~ 170 s in the X-ray 
light curve, which shows that a typical 'steep decay /plateau/normal decay' behaviour (Fig. 1). Swift's 
UltraViolet Optical Telescope (UVOT) did not detect a counterpart even though it started settled 
explosures only 77 s after the trigger 6 . A 2^m counterpart was detected with the United Kingdom 
Infra-Red Telescope (UKIRT, Hawaii) 20 min after the trigger 7 . Evidence that this burst occurred at 
high redshift, was given by the multi-band imager Gamma- Ray Burst Optical/Near-Infrared Detector 
(GROND, Chile) multiband imager (from g' band to K band), which indicated a photometric redshift 
of z = 8.0t°^ (ref. 7). 

We used the 3.6m Telescopio Nazionale Galileo (TNG, La Palma) with the Near Infrared Camera 
Spectrometer (NICS) and the Amici prism to obtain a low-resolution (R ~ 50) spectrum of GRB 090423 
~ 14 hrs after the trigger. NICS/Amici is an ideal instrument to detect spectral breaks in the continuum 
of faint objects because of its high efficiency and wide simultaneous spectral coverage (0.8-2.4 fim.). The 
spectrum (Fig. 2) reveals a clear break at a wavelength of 1.1 (ref. 8). We derive a spectroscopic 
redshift for the GRB of z = 8.1^1 (ref. 9; see Supplementary Information, section 3), interpreting the 
break as Lyman-a absorption in the intergalactic medium. No other significant absorption features were 
detected. This result is consistent, within the errors, with the measurement reported in ref. 7. 

At z ~ 8.1, GRB 090423 has a prompt-emission rest-frame duration of only Tg r f = 1.13 ±0.12 s in 
the redshifted 15-150 keV energy band, an isotropic equivalent energy E; LSO = 1.0 ±0.3 x 10 53 erg in the 
redshifted 8-1000 keV energy band 10 and a peak energy E Pj7 .f = 437 ± 55 keV. The short duration and 
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the high peak energy are consistent both with the distribution of long bursts, linked to massive stellar 

collapse, and with the population of short bursts, thought to arise from the merger of binary compact 

stars 11 ' 12 . Although the analysis of the spectral lag between the high- and low-energy channels in the 

BAT band is inconclusive about the classification of GRB 090423, the high Ei so argues in favor of a 

long GRB. The fact that GRB 090423 matches the E iso — E, p , r f correlation of long GRBs within 0.5(7 

further supports this classification 33 (Supplementary Fig. 2). 

The rest-frame 7— ray and X-ray light curve of GRB 090423 is remarkably akin to those of long 
GRBs at low, intermediate and high redshifts (Fig. 1), suggesting similar physics and interaction with 
the circumburst medium. The near-infrared light curve of GRB 090423 ~ 15h after the trigger shows 
a temporal decay with a power-law index of «o ~ 0.5, which is markedly different from the decay 
observed at X-ray energies during the same time interval, which has a power-law index of ax,2 ~ 1.3 
(Supplementary Fig. 3 and Supplementary Information, section 2). As for other lower-redshift GRBs, 
this behaviour is difficult to reconcile with standard afterglow models, although the sampling of the 
near-infrared light curve is too sparse for any firm conclusion to be drawn. 

The spectral energy distribution of near-infrared afterglow is well fitted by a power-law with an index 
of (3 = 0.4^4 and an equivalent interstella extinction of E(2? — V)< 0.15, assuming dust reddening 
consistent with the Small Magellanic Cloud 9 . On the other hand, the analysis of the XRT data in the 
time interval 3900s-21568s suggests the presence of intrinsic absorption (in excess of the Galactic value) 
with an equivalent hydrogen column density of Nh(z) = 6.8+53 x l^ 22 cm_2 (90% confidence level; 
Supplementary Information, section 1). The low value of the dust extinction coupled with a relatively 
high value of Nh suggests that GRB 090423 originates from a region with low dust content relative to 
those of low-z GRBs 14 , but one similar to that of the high-z GRB 050904, for which z = 6.3 (ref. 15). 
Because the absorbing medium must be thin from the point of view of "Thomson" scattering, the 
metallicity of the circumburst medium can be constrained to be > 4% of the solar value, Z Q . The 
implication is that previous supernova explosions have already enriched the host galaxy of GRB 090423 
to more than the critical metallicity, Z ~ 10~ 4 Z & (ref. 16) that prevents the formation of very massive 
stars (population III stars). Therefore the progenitor of GRB 090423 should belong to a second stellar 
generation. Its explosion injected fresh metals into the interstellar medium, further contributing to the 
enrichment of its host galaxy. Its existence empirically supports the cosmological models 17,18 in which 
stars and galaxies, already enriched by metals, are in place only ~ 600 million years after the Big 
Bang. Long GRBs are mostly associated with star forming dwarf galaxies, which are thought to be the 
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dominant population of galaxies in the early Universe 19 . The fact that GRB 090423 appears to have 
exploded in an environment similar to that of low-z GRB hosts 20 is in agreement with this. 

The occurrence of a GRB at z ~ 8 has important implications for the cosmic history of these 
objects 21 ' 22 ' 23,24 . In a first, simple approach, we can assume that GRBs trace the cosmic star formation 
history, given the well-known link of the long GRBs and the deaths of massive stars 25 , and that GRBs 
are well described by a universal luminosity function. However, under these assumptions the expected 
number of bursts at z > 8 with an observed photon peak flux larger than or equal to that of GRB 090423 
is extremely low: ~4x 10 -4 in ~ 4 yrs of Swift operation (Supplementary Fig. 6 and Supplementary 
Information, section 4). Hence, one or both of the above assumptions may be oversimplification 2 ' 26 . 
The detection of a very high-z burst such as GRB 090423 could be accomodated if the GRB luminosity 
function were shifted towards higher luminosity according to (1 + z) 5 with S ^ 1.5 or if the GRB 
formation rate were strongly enhanced in galaxies with Z ^ 0.2Zq. The requirement for evolution may 
be mitigate if we assume a very high star formation rate at z > 8. However, we note that the need for 
evolution is strongly supported by both the large number of Swift detections at z > 2.5 (ref. 24) and 
the number of bursts with peak luminosities in excess of 10 53 erg s _1 (ref. 26). A possible explanation is 
that high-redshift galaxies are characterized by a top-heavy (bottom-light) stellar initial mass function 
with a higher incidence of massive stars than in the local Universe 27 , providing an enhanced number of 
GRB progenitors. Such objects could be the main agents responsible for completing the reionization of 
the Universe 28 ' 19 ' 29 ' 30 . 
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Figure 1. Rest-frame 7-ray and X-ray light curves for bursts at different redshifts. BAT and 

XRT light curve of GRB 090423 (red data) in the source rest-frame. Errors on luminosity, Li SO , are at la level; horizontal 
bars refer to the integration time interval. The XRT 0.3-10 keV light-curve shows a prominent flare at a rest-frame time 
of t r f ~ 18 s (also detected by BAT), and a flat phase (with a power-law index of ax,i = 0.13 ± 0.11) followed by a rather 
typical decay with power-law index ax,2 = 1.3±0.1. We compare the light curve of GRB 090423 with those of seven GRBs 
in the redshift interval 0.8-6.3. The bursts are selected from among those showing a canonical three-phase behaviour (steep 
decay /plateau/normal decay) in the X-ray light curve and without a spectral break between BAT and XRT, allowing the 
spectral calibration of the BAT signal into the 0.3-10 keV energy band. The light curve of GRB 090423 does not show any 
distinguishing features relative to those of the lower-redshift bursts, suggesting that the physical mechanism that causes 
the GRB and its interaction with the circumburst medium are similar at z ~ 8.1 and at lower redshifts. 
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Figure 2. TNG spectrum of the NIR afterglow. Bottom Panel. Spectrum of GRB 090423 obtained using 
the Amici prism on the Telescopio Nazionale Galileo (TNG). The sharp break at wavelength A w 1.1 fim, which is due to 
the HI absorption in the intergalactic medium at the wavelength of the Lya line, implies that z = 8.1 jig \- The spectrum 
has been smoothed with a boxcar filter of width A = 25 pixels (where one pixel corresponds to ~ 0.006 /im at A = l.l^im). 
The absolute flux calibration was obtained by matching the almost simultaneous GROND photometric measurements 7 . 
The wavelength calibration was obtained from the TNG archive and adjusted to the wavelengths of the main atmospheric 
bands. The error bar corresponds to ±1<t uncertainty as measured on the smoothed spectrum. The confidence level of the 
Lyman-a break detection is £ 4<r. See also Supplementary Information, section 3. Top Panel. Plot of the trasmittance 
(the atmospheric transparency convolved with the instrumental response). The system has a significant sensitivity down 
to 0.9 /im, and no instrumental or atmospheric effect could explain the abrupt flux break observed in the spectrum of 
GRB 090423. 
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Supplementary Information 

This material presents technical details to support the discussion in the main paper. 
We discuss here the Swift data analysis, the analysis of the photometric data, the details 
of the analysis of the TNG spectroscopic data and the modelling of the GRB redshift 
distribution. 

1. Swift data analysis 

Swift-BAT triggered on GRB 090423 at 07:55:19.35 UT on 23 April 2009. BAT data were analysed 
using the HEASOFT software package (version 6.6.2) with the Swift Calibration Database (CALDB) 
version BAT(20090130). Background-subtracted light curves in different energy channels, energy spectra 
and corresponding response functions were derived from the BAT event file as processed with the BAT 
software tool batgrbproduct, by using the mask- weighting technique for the BAT refined position 4 , and 
by using standard and BAT-dedicated software tools. The mask-weighted light curve (Supplementary 
Figure 1) showed a couple of overlapping peaks starting at To — 2 s, peaking at To + 4 s, and ending 
at To + 15 s. The estimated duration, T90, was 10.3 ± 1.1 s for the mask-weighted light curve in the 
15-150 keV band. T90 is defined as the duration 90% of the total prompt 7-ray fluence in the observer 
frame (i.e. the interval from 5% to 95% of the total fluence) and is estimated using the battblocks 
software tool. Noteworthy is the light curve of the hardest channel, from 100 to 150 keV, showing a very 
weak signal as compared with those of the other energy channels. This reflects the spectral softness of 
this GRB, as shown also by the total energy spectrum. The latter was accumulated from -0.7 to 11.7 
s and is fitted with a cut-off power law, N(E) ~ T~ r exp[(2 — T)E/E p ], with the best-fit value for 
the peak energy, E p = 48^ keV, and a photon index T = 0.6^Q g. This value of E p is fully consistent 
with that determined by fitting the Fermi/Gamma-ray Burst Monitor spectrum 10 with the canonical 
"Band" function 31 . The corresponding total fluence in the 15-150 keV energy band is (5.9 ±0.4) x 10~ 7 
erg cm~ 2 . The 1-s peak photon flux measured from 3.5 s in the 15-150 keV band is 1.7 ± 0.2 ph s _1 
cm~ 2 . Uncertainties are given at 90% confidence. At z = 8.1, GRB 090423 is found to be consistent with 
the -Ep^f - Ei so correlation 13 within 0.5<r (Supplementary Figure 2). We note that, even considering the 
measured peak energy as obtained by a fit of the Fermi/GBM data with a cut-off power-law spectrum, 
i.e. E p = 82 ± 15 keV and thus E PtTf = 746 ± 137 keV (ref. 10 ), GRB 090423 would still be consistent 
within 2o" with the E p<r f — Ei so correlation. 

The BAT light curve shown in Supplementary Figure 3 is the mask-weighted curve extracted between 
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15 and 150 keV, binned so as to ensure S/N> 2 with a minimum binning time of 0.512 s. Extrapolation 
of the BAT flux down to the 0.3-10 keV band was performed by assuming the above spectral model. We 
note that the remarkable X-ray flare detected by XRT is seen in the BAT data as well. We note that 
assuming that the flare is still part of the prompt emission 32 , the total duration of the prompt phase in 
the source rest frame might be ~ 20s, similar to other long GRBs. 

The XRT observations began 73 s after the trigger: up to ~ 300 s the signal was dominated by 
a flare. As in many other GRBs, the light-curve then flattened to a shallow decay phase which could 
be well modelled by a power-law with index ax,i = 0.13 ± 0.11. At t ~ 4500 s the X-ray afterglow 
steepened to ax,2 = 1-3 ± 0.1 (errors at 68% confidence level). The flare was modelled by a standard 
profile 34 : this is characterised by a 1/e rise-time t r i se = 29.1 ± 3.6 s; 1/e decay-time tdecay = 65.5 ± 3.6 
s; 1/e width of At = 94.6 ± 7.3 s, while the asymmetry parameter is k = 0.38 ± 0.03. This implies a 
variability measure At/t pea k = 0.66 and a brightness contrast AFlux/Flux around 25. While the flare 
parameters are defined following ref. 34 , the reported uncertainties are worked out by using the entire 
covariance matrix. At the redshift of the burst, the flare has an energy Ei so = 3.6 x 10 51 erg in the 
redshifted 0.3-10 keV band of XRT, comparable to the energy released during the prompt emission of 
other GRBs. 

To evaluate the intrinsic column density absorbing the GRB 090423 spectrum, we extracted data 
in the 3900-21568 s time interval (observer frame). This interval was selected in order to avoid the 
bright X-ray flare whose variable spectrum might alter the fit and in order to have sufficient signal in 
the extraction region which we define as a count rate of more than 0.01 counts s -1 . The resulting 7984 
s exposure contains 680 counts in the range between 0.3-10 keV. The ancillary response file (arf) was 
created with the task xrtmkarf (within heasof t v.6.2.2) using the relevant exposure file and the latest 
v.ll reponse matrix function (rmf). The spectrum was binned to 20 counts per bin in order to assure 
a reasonable \ 2 statistic. 

We fit the X-ray spectrum with a composite absorption model consisting of a Galactic contribution 
and an intrinsic absorption fixed to z=8.1 using the tbabs model within the XSPEC (vl2.5.0aa) package. 
We left the Galactic value free to vary in the 2.9 — 3.2 x 10 20 cm" 2 range (based on the absorption maps 
by 35 and 36 ). The X-ray continuum was modeled with a power law, as is customary for the afterglow 
spectra of GRBs. The overall fit is good with a reduced Xred = 1-12 (28 degrees of freedom, corresponding 
to a null hypothesis probability of 30%). The resulting power law photon index is Tx = 1.97to'ig- 
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For the intrinsic column density, we get a value of Nh{z) = x 10 22 cm~ 2 (90% confidence 

level), among the highest of all Swift GRBs 37 . The results refer to a solar composition and metallicity. 
Assuming that the medium is not Thomson thick, a lower limit of the metallicity can be obtained by 
Nh{z) ~ (1/ 'gt)(Z I 'Zq)~ 1 , where cjt is the Thomson cross-section (e.g. ref. 38 ). We find Z > 0.043 Z & . 
A lower limit on the value of Nh(z) > 6 x 10 21 cm -2 is found at 95% confidence level corresponding to 
a lower limit on the metallicity of Z > 0.004 Z@. 

2. Analysis of the photometric data 

We analyzed all the available photometric data by using the /it zphrem code??, in order to de- 
termine the photometric redshift and spectral properties of the afterglow. Our code fits a model of 
functional form f u oc v~^t~ a ° , including dust extinction (by an SMC-type extinction law) as a free 
parameter. 

We find that the data are best fit within the time range 4.2 x 103 < t < 6.6 x 104 s) with a model 
characterised by a temporal decay with a power law index ao = 0.50 ± 0.05 (we quote hereafter 95% 
confidence intervals). The dust content is constrained to be E(B — V) < 0.15, and the photometric 
redshift is %>hot = 8.3 ± 0.3, consistent with the spectroscopic results. We caution that the rest-frame 
wavelength observed extends only out to 2500A, and only the three JHK filters do indeed measure any 
flux redwards of Lyman-a. That is the reason why the spectral index is only loosely constrained and its 
error is asymmetric ((3 = OA^'I), although its relatively blue color still enables us to put stringent limits 
to the possible dust content in the afterglow environment. Supplementary Figure 4 shows the projections 
of the (z,ao, (3, E(B — V)) four-dimensional confidence intervals on the different bidimensional planes. 

Extending the analysis to the whole available temporal window (2.5 x 10 2 < t < 1.4 x 10 6 s) renders 
impossible to find a good fit with a single temporal power-law, because of the different decay regimes 
that the afterglow goes through. 

3. Analysis of the TNG spectroscopic data 

We observed the afterglow of GRB090423 with the near-IR camera and spectrograph NICS 40 on 
the Italian 3.6m Telescopio Nazionale Galileo (TNG) at La Palma. We used the the lowest spectral 
resolution mode, offered by the Amici prism 41 . This prism provides a simultaneous spectral coverage 
over a wide wavelength range, between 0.8 and 2.4 /im, and has a high efficiency. It yields a constant 
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spectral resolution R « 50 over the whole wavelength range. These characteristics make the instrument 
especially well-suited for studying the spectral distribution of faint objects. 

We obtained 128 minutes of on-target spectroscopy. The afterglow was positioned in the slit us- 
ing as reference a nearby star approximately 30 arcseconds away (at J2000 coordinates 09:55:35.31, 
+18:09:03.9). We used a dithering mosaic of 8 cycles, each including two coadds of single 120s expo- 
sures, repeated 4 times. The mean time of our observations was Apr 23.98, approximately 15.5 hours 
after the burst detection. The 2-dimensional spectrum is shown in Supplementary Figure 5. 

Standard reduction tasks for NIR spectroscopy were performed independently by four different 
groups in our team, all of them reaching consistent results. Wavelength calibration was obtained by using 
a standard calibration table provided by the TNG and matching the deep telluric absorption bands. 
This method allows for wavelength calibrations better than 0.005 /im at 1.1 /im, and its contribution 
(Az = ±0.04) to the final error budget on redshift is negligible. 

Relative flux calibration was performed by using the observed spectral shape of the reference star. Its 
optical (SDSS) and near-IR (2MASS) colors are consistent with those of an M3-III star. The absolute 
calibration of the spectrum was obtained from the comparison with the simultaneous photometric 
measurements obtained by GROND (H=19.94 (Vega), ref. 7 ). We estimate the slit losses to be less than 
30%. 

The observed flux is compatible with zero below a wavelength of 1.1 fim, while a significant flux 
(> 99% confidence level) is measured redwards of this limit. Assuming that this is due to hydrogen 
absorption by a virtually completely thick Lyman-a forest, then the redshift at which the GRB occurred 
is z = 8.1+0:3- The quoted error includes the uncertainties on the wavelength calibration and on the 
estimate of the break position. This value makes GRB090423 the most distant object spectroscopically 
identified to date. By using a standard cosmology with f^A = 0.73, VLm = 0.27, Hq = 71 kms _1 Mpc _1 , 
we find that GRB 090423 was detected at a lookback time of greater than 13 Gyrs. 

We tentatively identified two absorption features at 1.3 and 2.2 fim. These would be consistent with 
blends of Si IV and Fe II at 1400A and 2400A, z = 8.1 rest-frame, respectively. The detection, however, 
has a low confidence level due to the low S/N of the spectrum. 

4. Modelling the GRB redshift distribution 

We compute the probability of detecting of GRB 090423 in three different scenarios for the formation 
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and cosmic evolution of long GRBs: (i) no evolution model, where GRBs follow the cosmic star formation 
and their luminosity function (LF) is constant in redshift; (ii) luminosity evolution model, where GRBs 
follow the cosmic star formation but the LF varies with redshift; (iii) density evolution model, where 
GRBs form preferentially in low-metallicity environments. In the first two cases, the GRB formation rate 
is simply proportional to the global cosmic star formation rate as computed by . For the luminosity 
evolution model, the typical burst luminosity is assumed to increase with redshift as (1 + z) s . Finally, 
for the density evolution case, the GRB formation rate is obtained by convolving the observed SFR 
with the fraction of galaxies at redshift z with metallicity below Z t h using the expression computed by 
43 . In this scenario, the LF is assumed to be constant. 

The computation works as follows (see also 21 > 22 > 23 > 24 > 26 > 4 3)_ The observed photon flux, P, in the 
energy band E m \ D < E < £ max , emitted by an isotropically radiating source at redshift z is 

4ndl(z) ' [l) 

where S(E) is the differential rest-frame photon luminosity of the source, and di{z) is the luminosity 
distance. To describe the typical burst spectrum we adopt the functional form proposed by 31 , i.e. a bro- 
ken power-law with a low-energy spectral index a, a high-energy spectral index (3, and a break energy 
Eb = (a — (3)E p /(2 + a), with a = — 1 and [3 = —2.25 (ref. 45 ). In order to broadly estimate the peak 
energy of the spectrum, E p , for a given isotropic-equivalent peak luminosity, L = Jiy^^ ^ ES(E)dE, 
we assumed the validity of the correlation between E p and L (ref. 46 ). 

Given a normalized GRB LF, 4>(L), the observed rate of bursts with P\ < P < P2 is 

dt 'Jo dz 4vr l + z Jl(p u z) 

where dV{z)/dz is the comoving volume element, AS1 S is the solid angle covered on the sky by the 
survey, and the factor (1 + z)' 1 accounts for cosmological time dilation, ^grb(^) is the comoving 
burst formation rate and the GRB LF is described by a power law with an exponential cut-off at low 
luminosities 47 , i.e. <p(L) cx (L/L cut )~^ exp(— L cut /L). 

For the three scenarios, we optimize the model free parameters (GRB formation efficiency, burst 
typical luminosity at z = and the power index £ of the LF) by fitting the differential number counts 
observed by BATSE (see ref. 24 ' 26 for a detailed description of the models and of the analysis). We 
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find that it is always possible to find a good agreement between models and data. Moreover, we can 
reproduce also the differential peak flux count distribution observed by Swift in the 15-150 keV band 
without changing the best fit parameters. On the basis of these results, we compute the probability to 
detect with Swift a GRB at z > 8 with photon flux P. The results are plotted in Supplementary Figure 
6 (top panels) together with the cumulative number of GRBs at z > 8 expected to be detected by Swift 
in one year of observations (bottom panels). From the plot it is clear that the no evolution model fails 
to account for the observation of GRB 090423, since only ~ 4 x 10~ 4 GRBs are expected to be detected 
at z > 8 in ~ 4 years of Swift observations. Evolutionary models (both in luminosity or in density) can 
easily account for the discovery of GRB 090423. We note that the results confirm the need for cosmic 
evolution in the GRB luminosity function and / or in the GRB density obtained by recent analysis of the 
whole Swift GRB dataset. Indeed, both the large number of z > 2.5 bursts 24 and the number of bright 
(i.e. with peak luminosity L > 10 53 erg s _1 ) bursts 26 strongly require the existence of evolution. 

Moreover, we want to stress here that our conclusions are conservative. First of all, many biases 
can hampered the detection of GRB at very high redshift. Indeed, a few very high-z bursts may be 
hidden among the large sample of Swift bursts that lack of an optical detection. Thus, the discovery of 
a single event at z > 8 in 4.5 yrs of Swift operation can be treated as a lower limit on the real number of 
high-z detection. Moreover, our choice of the GRB LF is also conservative, since the existence of large 
population of faint GRBs (i.e. for an LF with a more gentle decline or a rise in the faint end) would 
lead to a decrease of the expeceted number of GRBs at z > 8 strenghtening our conclusions. 
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Time since BAT trigger [s] 



Supplementary Fig. 1. BAT mask- weighted light curve. Four channels and combined 0.512 s mask- weighted 
light curve. The light curve of the 100-150 keV energy channel shows a weak signal, because of the soft spectrum; the 
corresponding integration time is 2.048 s. Errors are at la level. 



GRB 090423 at a redshift of z ~ 8.1 



17 




10 49 1 50 1 51 10 52 1 53 1 

E iso ( er §) 



Supplementary Fig. 2. Isotropic energy and peak energy correlation. Position of GRB 090423 in the E Pil{ 
- Ei BO plane based on Swift/BAT 4 and Fermi/GBM 10 (fitted with the Band function) results. The lines show the best-fit 
power-law and the ±2<j region of the correlation as derived by 13 . Also shown are the 70 GRBs included in the sample 
analyzed in that work (errors on individual bursts are at la level). Given that short GRBs do not follow the correlation 13 , 
this evidence supports the hypothesis that, despite its cosmological rest-frame duration of ~ 1.3 s, GRB 090423 belongs 
to the long GRB class. 
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Supplementary Fig. 3. Observed light curve. Light curve of GRB 090423 as observed by Swift/BAT (red 
crosses), Swift/XRT (blue plus) and in the NIR (cyan points). Errors on fluxes are at la level and horizontal bars refer 
to the integration time interval. The XRT 0.3-10 keV light-curve, starting at 73 s after the burst, shows a prominent flare 
at t ~ 170 s (also detected by BAT), and a fiat phase (ax,i = 0.13 ± 0.11) followed by a rather typical decay (starting 
at t = 4513 ± 491 s) with power-law index ax,2 = 1.3 ± 0.1. Available photometric data are plotted in the K band (AB 
magnitude) by transforming the fluxes, when the observations have been taken in a different filter, using a power law with 
/3 = 0.4, as estimated from the NIR spectral energy distribution. A small displacement in time for contemporary data 
in different bands is applied in order to increase the visibility. The NIR light curve is consistent with a plateau phase 
(t ~ 10 2 — 10 s s) followed by a decay with ao ~ 0.5 (t ~ 10 3 — 10 5 s). This decay phase is shallower than the X-ray 
decay in the same time interval. Triangles at t ~ 10 5 s report NIR upper limits as obtained by our second epoch TNG 
observation with the NICS camera in the Y and J band and by GROND in the JHK band. These limits are consistent 
with the temporal decay observed by XRT. 
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Supplementary Fig. 4. Multi-parameter analysis of the photometric data. Analysis of available photometric 
data for GRB090423 in the interval 70 min < t < 1100 min. The code fits a model function with temporal index ao and 
spectral index f3, dust extinction E(B — V), and redshift z. The different panels show the projection of the four-dimensional 
confidence intervals on the different two-dimensional planes of interest. The best-fit is marked by the black dot, with the 
red, cyan, and green contours defining respectively the 68%, 95%, and 99.5% confidence areas. The apparently bimodal 
distribution in the /3 direction is an artifact of the parameter space discretization. 




Supplementary Fig. 5. TNG 2-dimensional spectrum. The spectrum has been taken by ~ 14 hrs from the 
trigger. The spectrum of the nearby reference star is also shown. 
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Supplementary Fig. 6. Probability of the occurance of GRB 090423. Top panels: probability for a GRB 
with peak photon flux P to be detected by Swift at z > 8. Luminosity evolution models are shown in the left panel, where 
shaded area refers to a typical burst luminosity increasing as L cut cc (l + z) s with 8 = 1.5 — 3. Density evolution models are 
shown in the right panel, where shaded area refers to a metallicity threshold for GRB formation Z±h = 0.02 — 0.2 Zq (the 
lower bound refers to the higher Z t h)- In both panels, the dashed line shows the no evolution case. The red point marks 
the position of GRB 090423. We note that the point represents a lower limit on the number of detection at z > 8 since 
a few very high-z bursts may be hidden among those bursts that lack of an optical detection. Bottom panels: cumulative 
number of GRBs at z > 8 to be detected by Swift with photon flux larger than P in one year of Swift observations. 



